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ABSTRACT 



Aims. EX Lup is the prototype of EXors, a subclass of low-mass pre-main sequence stars whose episodic eruptions are attributed to 
temporarily increased accretion. In quiescence the optical and near-infrared properties of EX Lup cannot be distinguished from those 
of normal T Tau stars. Here we investigate whether it is the circumstellar disk structure which makes EX Lup an atypical Class II 
object. During outburst the disk might undergo structural changes. Our characterization of the quiescent disk is intended to serve as a 
reference to study the physical changes related to one of EX Lupi's strongest known eruptions in 2008 Jan-Sep. 
Methods. We searched the literature for photometric and spectroscopic observations including ground-based, IRAS, ISO and Spitzer 
data. After constructing the optical-infrared spectral energy distribution (SED), we compared it with the typical SEDs of other young 
stellar objects and modeled it using the Monte Carlo radiative transfer code RADMC. A mineralogical decomposition of the lO/jm 
silicate emission feature and also the description of the optical and near-infrared spectra were performed. 

Results. The SED is in general similar to that of a typical T Tauri star, though above 7 yum EX Lup emits higher flux. The quiescent 
phase data suggest low level variability in the optical-mid-infrared domain. Integrating the optical and infrared fluxes we derived 
a bolometric luminosity of 0.7 Lg. The lOyum silicate profile could be fitted by a mixture consisting of amorphous silicates, no 
crystalline silicates were found. A modestly flaring disk model with a total mass of 0.025 Mq and an outer radius of 150 AU was 
able to reproduce the observed SED. The derived inner radius of 0.2 AU is larger than the sublimation radius, and this inner gap sets 
EX Lup aside from typical T Tauri stars. 
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1. Introduction 

In 2008 January the amateur astronomer A. Jones announced 
that EXLupi, the prototype of the EXor class of pre- 
main sequence eruptive variables, had brightened dramatically 
(I Jones 2OO8I 1. The quiescent phase brightness of EX Lup, an 
MOV star in the Lupus 3 star forming region, is Va: 13 mag. 
During its unpredictable flare-ups its brightness may increase 
by 1-5 mag for a period of several months (Herbig 1977 1. In 
the last decades EX Lup produced a number of eruptions, the 
last one in 2002 ( |Herbig 1977HHerbig et al. 2001||Herbig 2007| l. 
During the present outburst, which lasted until 2008 September 
(AAVSO International Databas43) EX Lup reached a peak 
brightness of 8 mag in January, brighter than ever before. 

According to the current paradigm, eruptive phenomena in 
pre-main sequence stars (FU Ori and EX Lup-type outbursts) are 
caused by enhanced accretion onto the central star. During these 
phases of the star formation process intense build-up of the stel- 
lar mass takes place. (Hartmann& Kenyon 1996 1. In this pic- 
ture circumstellar matter must play a crucial role in the eruption 
mechanism. Surprisingly little is known about the circumstel- 
lar environment of the EXor prototype, EX Lup. While its in- 
frared excess emission detected by IRAS and ISO was attributed 
to a circumstellar disk ( Gras- Velazquez & Ray , 2005), no de- 
tailed analysis or modeling of the disk structure has been per- 
formed so far. Such an analysis could contribute to the clari- 
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fication of the eruption mechanism and also answer the long- 
standing open question of what distinguishes EXors from nor- 
mal TTau stars. jHerbig (2008)| found no optical spectroscopic 
features which would uniquely define EXors, and he also con- 
cluded that their location in the (J - H) vs. {H - Ks) color-color 
diagram coincides with the domain occupied by TTau stars. 

In this work we search the literature and construct a complete 
optical-infrared spectral energy distribution (SED) of EX Lup 
representative of the quiescent phase. The SED will be modeled 
using a radiative transfer code, and will be compared with typ- 
ical SEDs of T Tau stars, in order to reveal differences possibly 
defining the EXor class. During the recent extreme eruption EX 
Lupi was observed with a wide range of instruments. Our re- 
sults on the physical parameters of the system can be used as a 
reference for evaluating the changes related to this outburst. 



2. Optical / infrared data 

2.1. Optical and near-infrared pliotometry from the literature 

We searched the literature and collected available optical and 
near-infrared photometric observations of EX Lup obtained dur- 
ing its quiescent phase. We defined as quiescence periods when 
the source was fainter than 12.5 mag according to the visual 
estimates and V-band magnitudes in the AAVSO International 
Database. A visual inspection of the DSS and the 2MASS im- 
ages revealed no nebulosity surrounding the star, thus we con- 
cluded that measurements taken with different apertures can be 
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Table 1. Optical and infrared observations of EX Lup collected 
from the literature, obtained in its quiescent periods. References: 
(1) Bastian & Mundt (1979); (2) Mundt & Bastian (1980); 
(3) Herbig et al (1992); (4) The DENIS database (DENIS 
Consortium, 2005); (5) Glass & Penston (1974); (6) Hughes et 
al. (1994); (7) 2MASS All-Sky Catalog of Point Sources (Cutri 
et al. 2003) 
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safely compared. The query results are presented in Tab.[T]and 
plotted in Fig.[T] 

2.2. ESO 2.2/FEROS optical spectroscopy 

Three spectra of EX Lup in quiescent phase were taken in 2007 
July with FEROS at the 2.2 m MPG/ESO telescope in ESO La 
Silla, Chile. FEROS has a spectral resolut ion of 4 8000 and a 
wavelength coverage from 360 to 920 nm JKaufer et al. 1999). 
The data reduction was done using the online FEROS Data 
Reduction System implemented on the telescope, that produced 



39 one dimensional spectra of echelle orders and a merged spec- 
trum of the entire wavelength range. 

With a cross-correlation method we computed the projected 
rotational velocity and radial velocity of EX Lup in quiescence 
from the FEROS spectra. We cross-correlated the spectra with 
a numerical template mask jBaranne et al. 1996] l. The mean ra- 
dial velocity is -0.948 ± O.lkms"^ From the B-V color of 
EX Lup and by using a projected rotational velocity calibra- 
tion for FEROS ( |Setiawan et al. 2004 )| we derived a v sin / = 
4.4+2 km s"' . Sections of the spectrum showing the most promi- 
nent emission lines are plotted in Fig. 2. 

2.3. NTT/SOFI near-infrared spectroscopy 

We have found near-IR spectra in the ESO Data Archive. The 
spectra were obtained on 2001 May 4 using SOFI at the ESO 3.5 
m NTT telescope. The data were taken under the program 67. C- 
0221(A) (PI:D. FoUia). Two spectra were taken using the blue 
and red grisms (R ~ 1000), which cover the ranges 0.95 - 1.64 
and 1.53-2.52 yum, respectively. In both cases the total exposure 
time was 500 s, divided into 4 exposures obtained following an 
ABBA nodding pattern. The data were reduced using our dedi- 
cated routines developed under the IRAF environment. The data 
reduction included sky subtraction, wavelength calibration and 
frame combination. Spectra of the G3V star HIP 78466 were 
used to correct for telluric absorption. A modified version of the 
Xtellcor software jVacca et al. 2003l l was used in this step. The 
spectrum is presented in Fig. 3. 

2.4. Infrared Astronomical Satellite (IRAS) 

EX Lup is included in the IRAS Point Source Catalog with 
detections at 12, 25 and 60fj.m, and with an upper limit at 
lOOjum. The measurements took place during a quiescent phase 
of EX Lup. In order to refine these data and also try to obtain a 
lOOjum flux density value, we re-analyzed the IRAS raw scans 
using the IRAS Scan Processing and Integration tool (SCANPI) 
available at IPACQ. Depending on the wavelength, 6-8 scans 
crossed the object. The scans were re-sampled, averaged, and 
a peak was located within 30" of the position of EX Lup. Then a 
point-source template was fitted and fluxes were calculated. The 
flux errors were dominated by the uncertainty in estimating the 
background level. The fluxes were color-corrected in an iterative 
way by convolving the source SED and the IRAS filter profiles. 
The color-corrected fluxes and flux uncertainties can be seen in 
Tab.m 

2.5. Infrared Space Observatory (ISO) 

EX Lup was observed on 5 dates with ISOPHOT, the photometer 
on-board ISO (Le mke et al. 1996). Four measurements belong 
to a small monitoring programme which adopted similar instru- 
mental setup at each epoch (PI: T. Prusti). For such datasets an 
algorithm aimed to achieve high relative photometry accuracy 
was developed in our group dJuhasz et al. l2007). The raw data 
were corrected for instrumental effects following the scheme de- 
scribed in lJuhasz et al. 1 (2007). Absolute flux calibration was per- 
formed by using the on-board calibration source (FCS) for the 
3.6-12yL(m filters, and by assuming a default responsivity for the 
20-25 jjm range. At different single epochs 60, 100 and 200 yum 
high resolution maps were performed using the Astronomical 
Observing Template P32. These data were processed using a 

- http://scanpiops.ipac.caltech.edu:9000/applications/Scanpi/ 
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Fig. 1. Spectral energy distribution of EXLup, showing all data points from Tab. 1. and Tab. 2. and the smoothed spectra obtained 
with NTT/SOFl, ISOPHOT-S and Spitzer/IRS. The gray stripe marks the median SED of T Tauri objects from the Taurus-Auriga 
star-forming region. The 5-36 yum section of the median SED was constructed by Furlan et al. (20 06) based on Spitzer IRS data 
of 55 Class II objects with spectral types between K5 and M2, while the remaining part of our median SED was taken from 
ID'Alessio et al. I (1999) who computed their SED from 39 T Tauri stars. With dotted line we overplotted the stellar photosphere 
used in the models detailed in Sec. 4.3. 

Table 2. Color-corrected fluxes of EXLupi. For the ISOPHOT data we adopted a conservative formal uncertainty of 15%. 
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dedicated P32 Tool jTuffs & Gabriel 20031) : the calibration was 
performed using the related FCS measurements. The flux values 
were color-corrected by convolving the SED of EX Lup with the 
ISOPHOT filter profiles. The results are listed in Table 2, the 
quoted formal uncertainties correspond to a conservative esti- 
mate of 15%, however, the measurement error above 25 /vm is 
probably higher than this value. The fifth photometric dataset 
was published by Gras- Velazquez & Ray (2005), but these data 
will not be used in our further analysis due to their presumably 
higher uncertainties. 

At the four epochs of the monitoring programme also 
low-resolution spectrophotometry was performed with the 
ISOPHOT-S sub-instrument in the 2.5-1 1.6//m wavelength 
range. These observations were reduced following the process- 
ing scheme described by K ospal et al. (in prep., a)| which cor- 
rected for the slight off-center positioning of the source, sub- 
tracted the separate background spectra and adopted realistic 



error bars. The processed and calibrated spectra are shown in 
Fig.l. 



2.6. Spitzer Space Telescope 

IRAC. EX Lup was observed with IRAC at 3.6, 4.5, 5.8 and 
8.0//m on 2005 March 29 (PID: 3716, PI: G. Stringfellow), in 
the sub-array mode. At each of the nine dither positions 64 im- 
ages were obtained with an exposure time of 0.1 s. The BCD 
(S14.0) image cubes of the 64 frames were combined into two- 
dimensional images using the "irac-subcube-collapse" IDL rou- 
tine provided by the Spitzer Science Centre. Aperture photome- 
try was performed on the nine final images at each wavelength 
using a modified version of the IDLPHOT routines. The aper- 
ture radius was set to 3 pixels (3.66"), the sky background was 
computed in an annulus with an inner radius of 10 pixels and 
a width of 10 pixels (sigma clipping was used to reject outly- 
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Fig. 2. Details of the optical spectrum of EX Lupi showing significant or typical features. The interval shown in the second panel of 
the upper row is same as Herbig's (2007) Fig. 3. It is apparent that the broad component of the Fe II /i5018 line was absent in 2007. 
The third panel shows a typical region with narrow metallic lines. The O I A 7773 triplet, shown in the third panel of the lower row, 
was observed in absorption in several T Tauri stars (Hamann & Persson 1992). 



ing pixel values). The aperture correction was performed using 
the quoted values from the IRAC Data Handbook (V3.0, here- 
after IDH). Following the outline of the IDH (see Chapter 5) in 
the 3.6 fj.m band a pixel-phase correction was also applied to the 
measured flux densities. The flux density values measured in the 
nine dither positions (at each band) were averaged to get the final 
photometry. To calculate the final uncertainty we added quadrat- 
ically the measurement errors (obtained from the nine individual 
flux density values in each band) and an absolute calibration er- 
ror of 2% (IDH). 



IRS. EX Lup was observed with the Infrared Spectrograph 
(IRS, Houck et al. 2004) of the Spitzer Space Telescope on 
2004 Aug 30 (PID: 172, PI: N. Evans) and on 2005 March 
18 (PID:3716, PI: G. Stringfellow). On the first date the tar- 
get was measured using Short Low (5.2-14.5;um), Long Low 
(14.0-38 yum). Short High (9.9-19.5 A<m) and Long High (18.7- 
37.2/vm) modules. The integration time was 14 s for the low- 
resolution modules, while 30 s and 60 s were used for the Short 
High and Long High modules, respectively. At the second epoch 
only the Short Low, Short High and Long High modules were 
used and a PCRS peak-up was executed prior to the spectro- 
scopic observation to acquire the target in the spectrograph slit. 
The integration times was 14 s for the Short Low module with 4 
observing cycles for redundancy, while 120 s and 60 s were used 
for the Short High and Long High modules, respectively, both 
with 2 observing cycles. 

The spectra are based on the droopres and rsc prod- 
ucts processed through the S 15.3.0 version of the Spitzer data 
pipeline for the low- and high-resolution data, respectively. For 
the low-resolution spectra the background has been subtracted 
using associated pairs of imaged spectra from the two nodded 
position, also eliminating stray light contamination and anoma- 
lous dark currents. For the high-resolution spectra the back- 
ground has been removed by fitting a local continuum under- 
neath the source profile. Pixels flagged by the data pipeline 
as being "bad" were replaced with a value interpolated from 
an 8 pixel perimeter surrounding the flagged pixel. The spec- 
tra were extracted using a 6.0 pixel and 5.0 fixed-width aper- 



ture in the spatial dimension for the Short Low and the Long 
Low modules, respectively, while in case of the high-resolution 
modules the spectra were extracted by fitting the source profile 
with the known PSF in the spectral images. The low-level fring- 
ing at wavelengths >20//m was removed using the irsfringe 



package (Lahuis & Boogert 2003 1. The spectra are calibrated 



using a spectral response function derived from IRS spectra and 
MARCS stellar models for a suite of calibrators provided by the 
Spitzer Science Centre. To remove any effect of pointing offsets, 
we matched orders based on the point spread function of the 
IRS instrument, correcting for possible flux losses. The spectra 
obtained at the two diff'erent dates agree within an uncertainty of 
< 1 1%, and are shown in Fig.[T] 



3. Results 

3. 1 . Variability in fhe quiescent pliase 

The observations of EX Lup are sporadic, thus no complete SED 
could be constructed for any given epoch. We decided to merge 
all quiescent data regardless of their dates. In Fig. [T] we plotted 
all fluxes listed in Tabs.[T]and|2] At most wavelengths the scatter 
of data points is in the order of 25%, though at optical and a few 
mid-infrared wavelengths the difference between the lowest and 
highest values can be considerably larger The ISOPHOT-S and 
Spitzer spectra shown in Fig. [1] also exhibit differences at sim- 
ilar level. Although part of this scatter is related to instrumen- 
tal and calibration effects, the individual error bars are typically 
significantly smaller than the scatter of the data points, strongly 
suggesting the existence of an intrinsic variability. However, in 
case of the 60 and 100 yum data points the uncertainty of the mea- 
surements is much higher than at shorter wavelengths, thus here 
the difference between the fluxes is within the error bars. We 
consider the 25% peak-to-peak variation as an upper limit for 
the variability of EX Lup in quiescence. According to Fig. [1] 
the variability might be somewhat larger at optical wavelengths, 
but the low number of data points prevents us from claiming a 
wavelength-dependence. 
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Table 3. Observed wavelength, EW, and FWHM of the most 
prominent emission hnes in the optical spectrum of EX Lupi in 
quiescence* 
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23.0 
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*The full list of the emission lines identified in the spectrum are found 
in the electronic version of this Table. 



3.2. The optical-infrared SED 

The measured SED of EX Lup is presented in Fig.[T] The optical 
part is clearly dominated by the stellar photosphere. The optical 
and near-infrared colors, however slightly differ from the stan- 
dard colors of an MO V star. This fact was aheady mentioned 
by |Gras- Velazquez & Ray | (2005) who could not derive a pos- 
itive extinction value from the Eb-v and colors. Similarly 
Herbig et al. 200l| claimed it to be unknown as well. An infrared 
excess above the photosphere is detectable longwards of the K- 
band. The 3-8 fim range is smooth and devoid of any broad spec- 
tral features, indicating that EX Lup is neither deeply embedded 
to exhibit ices, nor hot enough for exciting PAHs. At 10//m a 
strong silicate emission appears, and a corresponding - though 
broader and shallower - silicate band can be seen around 20 fim. 
At longer wavelengths the continuum emission decreases fol- 
lowing a power law with a spectral index of about -4/3. Due to 
the lack of any sub-mm or mm measurement, the SED cannot be 
followed longwards of 200 fj.m. 



3.3. Spectroscopy 

3.3.1. The optical spectrum 

The high-resolution optical spectrum of EX Lupi, obtained on 
2007 July 30, shows a large number of emission lines in addi- 
tion to the photospheric absorption features characteristic of the 
young late-type star The emission lines are mostly symmetric, 
without noticeable P Cygni or inverse P Cygni absorption. The 
most prominent emission lines are the Balmer lines of the hy- 
drogen, the H and K lines and the infrared triplet of the ion- 
ized calcium, and the helium lines at 5875 and 6678 A. The 
observed wavelengths and equivalent widths in A, and the full 
widths at half maximum in kms"', determined by gaussian fit- 
ting, are listed in Table[3] In addition to these prominent features, 
more than 200 weak, narrow metallic lines could be identified 
using Moore's (1945) multiplet tables. The widths of the lines 
are around 10-20kms indicating that they originate from the 
active chromosphere of the star (Hamann & Persson 1992). All 



the emission lines identified in the spectrum are listed in the elec- 
tronic version of Tab. [3] Fig.|2]shows sections of the spectrum. 

The shape of the Ha line is nearly symmetric. Its equiva- 
lent width, W(Ha)=-35.9 A, is significantly larger than the up- 
per limit of the chromospheric Ha emission of MO type stars 
(~6 A, Barrado y Navascues & Martin 2003), indicating ac- 
tive accretion during the quiescent phase. The velocity width 
of the Ha line 10% above the continuum level is 362 km s"', 
larger than the lower limit of 270 km s"', set by White & Basri 
(2003) for accreting T Tauri stars. The empirical relationship 
between the 10% width of the Ha line and the accretion rate 
Mace, found by Natta et al. (2004), allowed us to derive Mace = 
4.2+2 3 ^ 10"'" Mo/yr Comparison of the Ha line with that pub- 
lished by Reipurth, Pedrosa & Lago (1996) shows that in 2007 
the Ha emission line of EX Lupi was more symmetric, narrower, 
and had somewhat larger equivalent width than in 1994 March, 
a few days after an outburst. 

3.3.2. Near-infrared spectrum 

The flux-calibrated infrared spectrum, observed on 2001 May 
4, is displayed in Fig. |3] The shape of the spectrum is similar 
to that of an unreddened M-type star without near-infrared ex- 
cess (Greene & Lada 1996). In the ZJ band the Paschen /3, y, 6, 
and the He I /1 1.083jum lines are seen in emission. Absorption 
features of the Na I at 2.21 //m, Ca I at 2.26 fim, and the rovibra- 
tional Av = +2 transitions of CO at 2.3-2.4- fim can be identified 
in the K-band part of the spectrum. The Bry line of the hydrogen 
is barely visible in emission, but its equivalent width cannot be 
measured. The following absorption features are also identified 
in our spectrum: at 1.199/im a very deep feature corresponding 
to S I; at 1.577 fim a feature associated to Mg I and at 1.673 fim 
another feature associated with Al I. 

The Pa/3 emission line is a useful accretion tracer: its lu- 
minosity correlates well with the accretion luminosity (e.g. 
MuzeroUe et al. 1998, Dahm 2008). We used the relationship 
established by Muzerolle et al. to derive the accretion lumi- 
nosity of EX Lup during quiescence from the measured flux 
F(Pay6)=1.7 X 10-'^ Wm-2. The resulting Lace = 3.7 x lO^^L© 
allows us to determine the mass accretion rate Mace, taking the 
mass and the radius of the star from Tab. 6. The accretion was 
assumed to proceed through a gaseous disk; for its inner radius 
we adopted a 5 R, (Gullbring et al. 1998). The result is 
Mace ~ 4.0 X 10"'° Mo/yr, in good agreement with Mace ob- 
tained from the velocity width of the Ha line. A flux calibration 
uncertainty of 20% and the uncertainty of the empirical relation- 
ship sets the lower and upper limit at 1.3 x 10"'° Mo/yr and 
1.1 X 10"'^ Mo/yr, respectively. 

4. Discussion 

4. 1 . Comparison witli young stellar objects 

Infrared SEDs of eruptive young stars (EXors and/or FUors) 
were examined by [Green et al. (2006)| [Quanz et aL] (2007) and 
Kospal et al. (in prep., b)' Some of these stars (eg. UZ Tau, VY 
Tau, DR Tau, FU Ori, Bran 76) resemble EX Lupi: their SEDs 
decrease towards longer wavelengths and show a silicate fea- 
ture in emission. The SEDs of the other group of FU Ori-type or 
EXor-like variables (eg. VI 057 Cyg, VI 647 Ori, PV Cep, OO 
Ser), exhibit flat or increasing SE Ds in the 20-100//m wave- 



length range (Ab raham et al. 20041 Kospal et al. 2007), which is 
clearly diff'erent from the shape of EX Lupi's spectral energy dis- 
tribution. [&een]erair(2006| and [Q^ suggested 
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A, ()a.m) 

Fig. 3. Near-infrared spectrum of EX Lupi, observed on 2001 
May 4. The region around the '"CO bandhead absorption is en- 
larged in the inset. 



that the diversity in the shape of the SEDs is related to the evo- 
lution of the system: younger objects, still embedded in a large 
envelope exhibit flat SEDs, while more evolved objects, emit- 
ting T Tauri-like SEDs have already lost their envelopes and 
have only circumstellar disks. According to this categorization 
we conclude that EX Lupi is relatively evolved among eruptive 
stars, and expect that its circumstellar environment consists only 
of a disk without an envelope. 

In order to compare EXLup with other - not eruptive - 
young stars, in Fig. [T] we overplotted a gray stripe marking the 
median SED of T Tauri objects from the Taurus-Auriga star- 
forming region (ID'Alessio et al. , 1999, Furlan et al. 2006). At 
shorter wavelengths the SED of EX Lupi follows the Taurus 
median. Longwards of approximately 7 yum, however, its abso- 
lute level becomes higher than the median by a factor of ~2.5. 
Nevertheless, even in this longer wavelength range its shape still 
resembles the Taurus slope. 

A further comparison of EXLupi with other pre-main 
sequence stars can be based on the mid-infrared spectrum. 
[Furlan et al. (2006)| presented a large sample of Spitzer IRS 
spectra of young stars and grouped them into several categories 
by measuring the slope of the SEDs and the strength of the 
10 and 20;um silicate features. A visual classification places 
EX Lup in their scheme in group A or in group B, defined by a 
relatively strong silicate feature and a flat or somewhat decreas- 
ing SED over 20 //m. [Furlan et al. (2006)| suggested that only 
limited dust-growth and settling has taken place in these groups. 
Nonetheless, we found that EX Lup is not a typical member of 
these groups, since shortwards of 8 yum the steep decrease of ex- 
cess emmission characteristic of their objects is missing in the 
case of EX Lupi. These quaUtative conclusions can be verified 
by calculating color indices for EX Lup. The strength of its sili- 
cate feature {Fiq - F com) IF com - 0.58 is an intermediate value 
between group A and group B objects. However, the n6-25 spec- 
tral index is typically negative for all objects in the sequence 
while positive (0.03) in the case of EXLupi, making it an 'out- 
lier' in the scheme. Also, the correlation between the «6_i3 and 
ni3_25 indices found by [Furlan et al. (2006)| does not seem to 
hold for EXLup. The colors of IRAS 04385+2550, a star ex- 
hibiting IRS spectrum very similar to that of EX Lup, were in- 
terpreted by [Furlan et al. (2006)[ as indication for the opening of 
an inner gap in the disk. This might also give an explanation 
in case of EX Lup for the lower 6 /im flux compared to that at 
longer wavelengths. Other sources from their sample that resem- 
ble EX Lup are UY Aur, CZ Tau and HP Tau. It is interesting to 



note that all four objects are binary systems with diff'erent sepa- 
rations. 

Finally we calculated the Thoi bolometric temperature and 
Lboi bolometric luminosity of EX Lup following the method of 
Chen et al. (1995). They studied the Taurus and Ophiuchus star 
forming regions, and analyzed the distribution of young stars of 
different evolutionary stages in the Lboi vs. Tboi diagram. In a 
later paper (Chen et al. 1997) they repeated the study for the 
objects of the Lupus cloud as well, including EX Lup. They 
found Tboi - 2365 K and Lhoi - 0.7 Lq based on IRAS data. 
Using our newly constructed SED we recalculated these two val- 
ues. Assuming Ay - Omag, our result is Tboi - 1982 K and 
Lboi - 0.73 Lq. According to these parameters, EX Lup seems to 
be a typical classical T Tauri star. The location of EX Lup on the 
Lboi vs. Thoi diagram implies that it is a Class 11 object, and its 
age is 3.8+^;^ X 10* yr. 

4.2. Mineralogy 

Based on the previous results that EX Lup is similar to classical 
T Tauri stars, in the following we assume that its circumstellar 
matter is distributed in a disk. In order to derive the dust com- 
position in the surface layer of the disk we fitted the Spitzer 
IRS spectrum using the Two-Layer Temperature Distribution 
(TLTD) method (Juhasz et al. 2009). This method assumes 
multi-component continuum (star, inner rim, disk midplane) and 
a distribution of temperatures to fit the source function. The sili- 
cate emission is expected to arise from the disk atmosphere. The 
observed spectrum is computed as 
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Here, d is the distance, By(T) is the Planck-function and k/j is 
the mass absorption coefficient of dust species / and grain size 
j. The subscripts 'rim', 'atm' and 'mid' denote the quantities 
of the inner rim, the disk atmosphere and the disk mid-plane, 
respectively. The dust species used are listed in Tab. 4. 

In Fig. [4] we present the fit to the 5-17jum region of the 
Spitzer IRS spectrum. The derived dust mass fractions are given 
in Tab. 5. The spectral decomposition shows that the main con- 
tributors to the optically thin lOjum silicate emission complex 
are the amorphous silicates (98.4 ± 1.0% in terms of mass). The 
mass fraction of crystalline silicates is below 2 %, which agrees 
with the value for the diffuse ISM (Kemper et al. 2005 1. It is in- 
teresting to note, that the mass-weighted average grain size of the 
amorphous silicates (0.57 jum) in our fit is somewhat larger than 
typically found in the ISM. At wavelengths longer than 17jum 
the IRS spectrum of EX Lup shows only the 18 /urn feature of 
the amorphous silicates, no crystalline feature can be seen. 

The fitted values for Trim.max, T'atm.max, T'mid.max are 1202±4K, 
1032 ± 90K and 846 + 4K, respectively. The fitted values for 
qnm, qmm, ?mid are -2.02 ± 0.01, -0.44 ± 0.02 and -0.42 ± 01, 
respectively. Although the;^'^ of the fit is 14.5, which is far higher 
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Fig. 4. Fit of the observed spectrum. The spectral decomposi- 
tion shows that the mid-infrared spectrum of EX Lup can be re- 
produced by a mixture of amorphous siHcates with oUvine and 
pyroxene stoichiometry. The contribution of crystalline silicates 
derived from the spectral decomposition is below 2%. 



than unity, expected for a good fit, the average deviation from 
the observed spectrum is about 1 %. The reason for the high;i'^ 
can be found in the realtively high signal-to-noise ratio of the 
Spitzer IRS spectrum (>300 in the fitted wavelength range) and 
in the deficiencies of the applied dust model (optical constants 
and grain shape model). The latter is probably responsible for 
the differences between the 2D RT model and the observed IRS 
spectrum between 14 and 21 jjm. 

Given the fact that EX Lup is a young eruptive star, one 
would expect to observe an increased value of crystallinity in 
the mid-infrared features, compared to 'normal' TTauri stars. 
Since the only requirement of the crystallization is the high tem- 
perature, the enhanced irradiation luminosity and viscous heat- 
ing during the outbursts should lead to rapid crystallization in 
the disk. lQuanz et al.| (2007), in accordance with our results, re- 
ported the lack of crystalline emission and the presence of emis- 
sion from larger grains (a > 0.1 /im) in the mid-infrared spectra 
of FU Orionis objects. They hypothesized that the reason for the 
lack of crystals can be twofold. One can be the replenishment of 
the dust content of the disk atmosphere (where the mid-infrared 
features originate) by pristine dust from an infalling envelope of 
the FU Ori object. Another possible explanation is strong verti- 
cal mixing in the disk which transports the crystals into deeper 
layers of the disk where they cannot be detected any more by 
mid-infrared spectroscopy. Out of these two speculative scenar- 
ios the latter is favorable in case of EX Lup since, to our knowl- 
edge, there is no infalling envelope around the source. 

4.3. Modeling 

In this section we perform a detailed modeling of the EX Lup 
system in order to derive the geometry of its circumstellar en- 
vironment. We fit the data points presented in Fig. 1, except the 
mid-infrared domain where, due to the intrinsic variability and 
the different quality of the data, we considered the Spitzer mea- 
surements only: the IRAC and IRS observations form a quasi- 
simultaneous high quality data set covering the 3.6 - 38 fim 
wavelength range. Moreover, these measurements are close in 



Table 4. Overview of dust species used in fitting the 5-17//m 
spectrum. For each component we specify its lattice struc- 
ture, chemical composition, shape and reference to the labora- 
tory measurements of the optical constants. For the homoge- 
neous spheres we used Mie theory to calculate the mass absorp- 
tion/scattering coefficients. For the inhomogeneous spheres, we 
used the distribution of hollow spheres dMin et al. 200 5), to sim- 
ulate grain shape deviating from perfect symmetry. References: 
(1) l^rsc hner et al. (1995), (2) [Servoin & Piriou (1973)| (3) 
JageretaL(1998), (4) |Henning & Mutschke (1997)| 



Species 


State 


Chemical 
formula 


Shape 


Ref 


Amorphous silicate 


A 


MgFeSi04 


Homogeneous 


(1) 


(olivine stoichiometry) 






sphere 




Amorphous silicate 


A 


MgFeSijOfi 


Homogeneous 


(1) 


(pyroxene stoichiometry) 






sphere 




Forsterite 


C 


Mg2Si04 


Hollow sphere 


(2) 


Chno Enstatite 


C 


MgSiOi 


Hollow sphere 


(3) 


Silica 


A 


SiO, 


Hollow sphere 


(4) 



Table 5. Fitted dust composition to the Spitzer IRS spectrum. 
We only show mass fractions that are larger than 0. 1 %. 

Dust species Mass fraction [%] 

O.lyum 1.5^m 6.0 jim 



Am. Silicate 


65.6 ± 0.5 




(Olivine type) 






Am. Silicate 




32.80 ± 0.5 


(Pyroxene type) 






Forsterite 


0.7 ± 0.08 




Enstatite 




0.5 ± 0.3 0.4 ± 0.4 


Silica 







time to the recent 2008 eruption and thus document the pre- 
outburst conditions of the system. 

For the stellar parameters (radius effective temper- 
ature Tt and mass M*) we adopted the values used by 
(2005). 



Gras- Velazquez & Ray (2005). These parameters were fixed 
during the modeling (Tab. 6 lists all input parameters marking 
differently the fixed and the optimized ones). For the distance 
of the star we adopted 155 pc, the mean distance of the Lupus 
Complex as measured by the Hipparcos dLombardi et al. 20081 1. 
Note that the unknown location of EX Lup within the complex 
introduces an additional uncertainty of about 30 pc. For the value 
of extinction we assumed Ay - mag (Sect. 3.2). 

We used the Monte Carlo radiative transfer code RADMC 
(IDullemond & Dominik 2004] l combined with RAYTRAC^ 
The circumstellar environment is supposed to be axially sym- 
metric, so we use a two-dimensional geometry in polar coordi- 
nates (r,ff}. As mentioned in Sects. 2.1 and 4.1 we had no indi- 
cation for the presence of an envelope around EX Lupi, thus our 
system consists only of a central star and a dusty disk. For the 
structure of the disk we assume the density profile as follows: 



Pdiskir, z) = 



nr) 1 

exp < — 

/z(r)V2^ I 2 



h{r) 



for details see: |http://www.mpia.de/homes/dulIemon/radtrans/radmc/| 
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where r and z are the radial and vertical coordinates respectively, 
2(r) = '^disk(r^) ^ is the surface density. The 'disk' subscript 

\ ' disk / 

denotes values of the given parameter at the outer radius. The 
scale height h(r) increases with radius as: 

r rdisk \rdiskj 

where Ofi is the flaring index. A schematic picture of the used 
geometry is presented in Fig.|5] 

Dust + gas disk Star 

Fig. 5. Schematic picture of the geometrical structure of the 
model. The figure is not to scale. 

The temperature distribution is determined by the heating 
sources: the central star, for which we used a Kurucz model at- 
mosphere (Castel ir& Kurucz 20031 1 and the heated dust grains 
emitting blackbody radiation. We used a passive disk and did 
not consider accretion, due to the low accretion rate in quies- 
cence (Sect. 3.3). The accretion luminosity is < 1% of the stellar 
luminosity, thus its contribution is negligible compared to direct 
irradiation coming from the central star. After the calculation of 
the temperature distribution, the SED of the system is produced 
at an inclination angle § with the ray-tracer. 

We used the dust composition derived from the fitting of 
the Spitzer IRS spectrum (Sect. 4.2), but excluded dust species 
whose contiburtion to the total mass is < 1%. This way, the dust 
model used in RADMC contained only amorphous silicates of 
olivine and of pyroxene types with a mass ratio of 2:1. The size 
of the dust grains were 0. 1 yum and 1 .5 //m for olivine and pyrox- 
ene stoichiometry, respectively. Besides these, we added 20 % 
amorphous carbon with a grain size of 0. 1 /im. The mass ab- 
sorption coeflicients of amorphous carbon were calculated using 
Mie-theory from the optical constants of lPreibisch et al.l (1993). 

We could fit reasonably most of the measured data points 
by adopting the geometry described above. In the best model of 
this type we had to move the inner radius of the dusty disk out to 
0.5 AU, significantly exceeding the dust sublimation radius (at 
T=1500K) of less than 0.1 AU. However, in the 3-8 /^m wave- 
length range this model underestimated the measured points 
(Fig lSl dotted line). Using smaller values for the inner radius, 
we could improve the fit of the 3-5 //m range, but then it was not 
possible to reproduce the measured fluxes at longer wavelengths 
(Fig lSl dashed line). Reducing the inner radius below 0.2 AU 
the model fluxes even in the near-infrared wavelength range be- 
came too high. In order to solve this problem, we introduced a 
rounded inner rim to the disk (see inset in Fig.|5]l, thus instead of 
having a sharp inner edge with a high wall we decreased the disk 
height to - 0.05 at the beginningof the disk and at the same 

' disk 

time moved the inner radius inward. Behind 0.6 AU we left the 
structure of the disk unchanged. With this modification we could 
obtain a fit that is in good agreement with all infrared observa- 
tions. Nevertheless, the inner radius of 0.2 AU used in the best 
model is still beyond the sublimation radius. 



Table 6. Parameters used in the best-fit model. Those in italics 
were adopted from the literature and kept fixed during the mod- 
eling. 



Parameters Fitted value 

System Parameters 

Distance (d) 155 pc 

Inclination {&) 20 ° 

Stellar Parameters 

Temperature (Ts,ar) 3800 K 

Mass (Msiar) 0.6 Mo 

Radius (Rstar) 1-6 R© 

Visual extinction (Ay) mag 
Circumstellar Disk Parameters 

Inner radius of dusty disk (r„, ,;,jt) 0.2 AU 

Outer radius of dusty disk (rj, jj.) 150 AU 

Scale height (^) 0.12 

' disk 

Flaring index (cfyv) 0.09 

Exponent of radial density profile (p) -1.0 

Total mass (M) 0.025 Mq 



















> 


; i 
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< 1 1 1 





1 10 100 
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Fig. 6. Spectral energy distribution of EX Lup. The solid line 
shows our best fit model with the rounded inner wall, while the 
dashed and dotted lines coiTespond to the best models using a 
sharp inner edge. The plotted error bars mark the range of ob- 
served fluxes at different epochs, taking into account also the in- 
dividual measurement uncertainties (filled dots correspond to the 
middle of the range). Error bars on the IRAC points are smaller 
then the size of the symbols. 



Our best-fit model is presented in Fig. |6] with a solid line, the 
coiTesponding parameters are listed in Tab.|6] In the literature we 
could not find any constraints for the inclination of the system 
and its value is not well defined by the model either. The disk is 
definitely seen closer to face-on then edge-on, and we obtained 
our best fit using a value of 20°, though fits for inclinations be- 
tween 0° (face on) and 40° are very similar. The outer radius 
of the circumstellar disk is typical for young systems, while the 
disk is massive among T Tauri disks. For the value of the flaring 
index a/i = 0.09 gave the best result, which means that the disk 
flares only very modestly. The scale height is /i = 12 AU at the 
outer radius of the disk, which is slightly higher than the value 
typical for T Tauri disks. The exponent of the radial surface den- 
sity profile is p = -1.0. 
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5. Summary and conclusions 

We characterized the quiescent disk of EXLupi and investi- 
gated whether it is the structure of the circumstellar environment 
which makes it an atypical, eruptive young stellar object. Our 
main findings are the following: 

- During quiescent phase there is indication for an intrin- 
sic variability of less than 25% in the optical-mid-infrared 
wavelength regime. 

- Our new spectra is consistent with the classification of 
EX Lup as an M type star Based on the Ha and Pa^ spec- 
tral lines, we derived a very low quiescent accretion rate of 
~4x lO-'^Mo/yr 

- In general the shape of the SED is similar to those of typi- 
cal T Tau stars, but above 7 fim EX Lup is brighter than the 
Taurus median by a factor of ~2.5. The relative flux contri- 
bution from shorter and longer wavelengths is a parameter 
which may distinguish EX Lup from the majority of classi- 
cal TTau stars. 

- The lOyum silicate feature of EX Lup can be well reproduced 
by amorphous silicates with olivine and pyroxene stoichiom- 
etry, no crystalline silicates were found. 

- A modestly flaring disk model with a total mass of 0.025 Mq 
and with inner and outer radii of 0.2 and 150 AU, respec- 
tively, is able to reproduce the observed SED. The radius of 
the inner hole is larger than the dust sublimation radius. 

The existence of this dust free inner hole points to a clear- 
ing mechanism, for which several explanations could be in- 
voked: (l) |Eisner et al. (2007)| claimed that low mass young stars 
with small accretion luminosities tend to have inner disk radii 
larger than the sublimation radius, probably due to efl'ects of 
the magnetic field. Using the relationship in Eis ner et al. (2007)| 
and taking 0.2 AU as the magnetospheric inner radius and an 
accretion rate of 4 x 10"'"Mo/yr EXLupi should have a mag- 
netic field strength of 2.3 kG, which is typical for T Tauri stars 
dJohns-KrullIOOTl l. 

(2) Binarity might also be responsible for clearing up regions of 
the disk. In case of EX Lup binarity was studied by several au- 
thors. |Ghe£etair(T997) used high angular resolution techniques 
in order to find wide components, and detected no companion 
of EX Lup between 150-1 800 AU separation. Bailey (1998) de- 
tected no companion between 1 - 10 AU. Melo (2003) claimed 
that EX Lup is not a spectroscopic binary, and similarly Herbig 
(2007) concluded the same based on spectroscopic data mea- 
sured by the Keck telescope. |Guenther et al. (2007)| performed 
an eight year long radial velocity monitoring program, but they 
could not detect binarity in case of EX Lupi either though they 
only had 3 spectra. Nevertheless, we note that if the inclination 
of EX Lupi is indeed close to a face-on geometry it makes detec- 
tion of binarity difficult. 

(3) Inner gaps are characteristic of disks in transitional phase be- 
tween the Class II and Class III stages. These transitional disks, 
however - unlike EX Lup - exhibit no excess over the photo- 
sphere in the near-infrared wavelength range. Furlan et al. (2006) 
suggested that IRAS 04385+2550, an object whose SED is sim- 
ilar to that of EX Lup, could be in a state preceding the transi- 
tion. Adopting this explanation, EX Lup might also be in a pre- 
transitional state. 

(4) Alternatively with such a low accretion rate photoevapora- 
tion by EUV radiation may also contribute to the clearing of this 
innermost region t iGorti & HoUenbach 20091 ). 

(5) It cannot be excluded that the large inner hole is a 



pheanomenon in connection with the eruptive behaviour, al- 
though details of this connection are unclear yet. 

Our detailed modeling of the quiescent disk structure of 
EX Lup (Sect. 4.2) may be a good basis for studying the physical 
changes related to the 2008 eruption. Assuming that the outburst 
is due to temporarily increased accretion, one could scale up the 
accretion rate in our quiescent model, and check whether this 
new model would reproduce the outburst SED. If this strategy 
fails to provide a sufficiently good fit to the outburst data, it may 
be a hint for geometrical restructuring of the circumstellar envi- 
ronment during the outburst. 

In the present work one of our main goals was to identify 
atypical features in the circumstellar structure of EX Lup which 
may explain its eruptive nature. The inner disk hole revealed by 
our modeling is an unexpected result in this sense and compari- 
son of EX Lup with other EXors at infrared wavelengths would 
be important. It may answer the question whether an inner gap 
in the dusty disk is characteristic of the EXor phenomenon thus 
the hole is connected to the eruption mechanism, and we could 
learn to which extent is EX Lup a good representative of eruptive 
stars. 

Acknowledgements. The work was supported by the grants OTKA T 49082 and 
OTKA K 62304 of the Hungarian Scientific Research Fund. The authors are 
grateful for C. P. Dullemond for kindly providing RAYTRACE and also for Jeroen 
Bouwman for providing his routines for the Spitzer IRS data reduction. The 
research of A. K. is supported by the Nederlands Organization for Scientific 
Research. They also thank for observational data from the AAVSO International 
Database contiibuted by observers worldwide. 



References 

Abraham, P., Kospal, A., Csizmadia, S., Moor, A., Kun, M., & Stringfellow, G. 

2004, A&A, 419, L39 
Bailey, I. 1998, MNRAS, 301, 161 
Baranne, A., et al. 1996, A&AS, 1 19, 373 
Barrado y Navascues, D., & Martm, E. L. 2003, AJ, 126, 2997 
Bastian, U., & Mundt, R. 1979, A&AS, 36, 57 

Castelli, F., & Kurucz, R. L. 2003, ModelHng of Stellar Atmospheres, 210, 20P 
Chen, H., Myers, R C, Ladd, E. F, & Wood, D. O. S. 1995, ApJ, 445, 377 
Chen, H., Grenfell, T. G., Myers, R C, & Hughes, J. D. 1997, ApJ, 478, 295 
Cutri, R. M., et al. 2003, The IRSA 2MASS All-Sky Point Source Catalog, 

NASA/IPAC Infrared Science Archive 
D'Alessio, P, Calvet, N., Hartmann, L., Lizano, S., & Canto, J. 1999, ApJ, 527, 

893 

Dahm, S. E. 2008, AJ, 136, 521 

The Denis Consortium 2005, VizieR Online Data Catalog, 1, 2002 
Dorschner, J., Begemann, B., Henning, Th., Jager, C, Mutschke, H. 1995, A&A, 
300, 503 

Dullemond, C. P, & Turolla, R. 2000, A&A, 360, 1187 
Dullemond, C. P, & Dominik, C. 2004, A&A, 417, 159 

Eisner, J. A., Hillenbrand, L. A., White, R. J., Bloom, J. S., Akeson, R. L., & 

Blake, C. H. 2007, ApJ, 669, 1072 
Furlan, E., et al. 2006, ApJS, 165, 568 

Ghez, A. M., McCarthy, D. W., Patience, J. L., & Beck, T. L. 1997, ApJ, 481, 
378 

Glass, I. S., & Penston, M. V. 1974, MNRAS, 167, 237 
Gorti, U., & Hollenbach, D. 2009, ApJ, 690, 1539 
Gras-Velazquez, A., & Ray, T. P 2005, A&A, 443, 541 

Green, J. D., Hartmann, L., Calvet, N., Watson, D. M., Ibrahimov, M., Furlan, 

E., Sargent, B., & FoiTest, W. J. 2006, ApJ, 648, 1099 
Greene, T. P, & Lada, C. J. 1996, AJ, 112, 2184 

Guenther, E. W., Esposito, M., Mundt, R., Covino, E., Alcala, J. M., Cusano, F, 

& Stecklum, B. 2007, A&A, 467, 1 147 
GuUbring, E., Hartmann, L., Briceno, C, & Calvet, N. 1998, ApJ, 492, 323 
Hamann, F, & Persson, S. E. 1992, ApJS, 82, 285 
Hartigan, P, & Kenyon, S. J. 2003, ApJ, 583, 334 
Hartmann, L., & Kenyon, S. J. 1996, ARA&A, 34, 207 
Henning, Th. & Mutschke, H. 1997, A&A, 327, 743 
Herbig, G. H. 1977, ApJ, 217, 693 

Herbig, G. H., Suntzeff, N., & Blanco, B. 1992, Infomiation Bulletin on Variable 
Stars, 3755, 1 



10 



N. Sipos et al.: EXLupi in quiescence 



Herbig, G. H., Aspin, C, GUmore, A. C, Imhoff, C. L., & Jones, A. F. 2001, 

PASP, 113, 1547 
Herbig, G. H. 2007, AJ, 133, 2679 
Herbig, G. H. 2008, AJ, 135, 637 

Hughes, J., Hartigan, P., Krautter, J., & Kelemen, J. 1994, AJ, 108, 1071 
Jager, C, Molster, F. J., Dorschner, J., Kenning, Th, Mutschke, H., Waters, 

L.B.FM., 1998, A&A, 339, 904 
Johns-Krull, C. M. 2007, ApJ, 664, 975 

Jones, A. F. A. L. 2008, Central Bureau Electronic Telegrams, 1217, 1 
Juhasz, A., Prusti, T., Abraham, P, & Dullemond, C. P 2007, MNRAS, 374, 
1242 

Juhasz, A., Henning, Th., Bouwman, J., Dullemond, CP, Pascucci, I., Apai, D. 

2009, ApJ, accepted 
Kaufer, A., Stahl, O., Tubbesing, S., N0rregaard, P., Avila, G., Francois, P., 

Pasquini, L., & Pizzella, A. 1999, The Messenger, 95, 8 
Kemper, F, Vriend, W., Tielens, A. 2005, ApJ, 633, 534 

K6spal, A., Abraham, P., Prusti, T., Acosta-Pulido, J., Hony, S., Mo(5r, A., & 

Siebenmorgen, R. 2007, A&A, 470, 211 
K6spal, A., Nemeth, P, Abraham, P., Kun, M., Henden, A., & Jones, A. F. 2008, 

Information Bulletin on Variable Stars, 5819, 1 
K6spal A., Abraham, P., Kun, M., Leinert, C, Henning, Th., & Acosta-Pulido, 

J. 2009a in prep. 
Kospal, A. et al. 2009b, in prep. 

Lahuis, F., & Boogert, A. 2003, SFChem 2002: Chemistry as a Diagnostic of 

Star Formation, 335 
Lemke, D., et al. 1996, A&A, 315, L64 
Lombardi, M., Lada, C. J., & Alves, J. 2008, A&A, 480, 785 
Melo, C. H. F 2003, A&A, 410, 269 

Min, M., Hovenier, J., & de Koter, A. 2005, A&A, 432, 909 
Moore, C. E. 1945, Contrib. Princeton Univ. Obs., No. 20 
Mundt, R., & Bastian, U. 1980, A&AS, 39, 245 
Muzerolle, J., Hartmann, L., & Calvet, N. 1998, AJ, 116, 2965 
Natta, A., Testi, L., Muzerolle, J., Randich, S., Comer6n, R, & Persi, P. 2004, 
A&A, 424, 603 

Preibisch, T., Ossenkopf, V., Yorke, H., Henning, Th. 1993, A&A, 279, 577 
Quanz, S. P., Henning, Th., Bouwman, J., van Boekel, R., Juhasz, A., Linz, H., 

Pontoppidan, K. M., Lahuis, F, 2007, ApJ, 668, 359 
Reipurth, B., Pedrosa, A., & Lago, M. T. V. T. 1996, A&AS, 120, 229 
Servoin, J. L. & Piriou, B. 1973, Phys. Stat. Sol., 55, 677 

Setiawan, J., Pasquini, L., da Silva, L., Hatzes, A. P., von der Liihe, O., Girardi, 

L., de Medeiros, J. R., & Guenther, E. 2004, A&A, 421, 241 
Tuffs, R. J., & Gabriel, C. 2003, A&A, 410, 1075 
Vacca, W. D., Gushing, M. C., & Rayner, J. T. 2003, PASP, 115, 389 
White, R. J., & Ghez, A. M. 2001, ApJ, 556, 265 
White, R. J., & Basri, G. 2003, ApJ, 582, 1 109 



